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F.G.S.L. Brandao, 1 ' 2 M.J. Hartmann, 1 - 2 and M.B. Plenio 1 ' 2 

1 QOLS, Blackett Laboratory, Imperial College London, London SW7 2AZ, UK 
2 IMS, Imperial College London, 53 Prince's Gate Exhib. Rd, London SW7 2PG, UK 

(Dated: February 1, 2008) 

We propose a new method to produce self- and cross-Kerr photonic nonlinearities, using light- 
induced Stark shifts due to the interaction of a cavity mode with atoms. The proposed experimental 
set-up is considerably simpler than in previous approaches, while the strength of the nonlinearity 
obtained with a single atom is the same as in the setting based on electromagnetically induced 
transparency. Furthermore our scheme can be applied to engineer effective photonic nonlinear 
interactions whose strength increases with the number of atoms coupled to the cavity mode, leading 
to photon-photon interactions several orders of magnitude larger than previously considered possible. 
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Quantum properties of light, such as photon anti- 
bunching [l[ and photonic entanglement 0], can only 
be produced by nonlinear interactions between photons. 
Strong nonlinearities are also important for quantum 
information processing, with applications ranging from 
quantum nondemolition measurements 0] to quantum 
memories for light [4] and optical quantum computing 
architectures [5(. However, photon-photon interactions 
are usually extremely weak, and several orders of mag- 
nitude smaller than those needed in the above applica- 
tions. A possible route towards larger nonlinearities is 
the use of coherent interaction between light and matter 
in high finesse QED cavities [6fl . The goal here is to pro- 
duce large nonlinearities with negligible losses, something 
which cannot be accomplished by merely tuning close to 
resonance atom-light interactions. To the best of our 
knowledge, the setting generating the largest nonlinear 
interactions in this context, proposed by Imamoglu and 
co-workers [l, Hi, [HJ , uses quantum interference effects re- 
lated to electromagnetically induced transparency (EIT) 
Q to produce giant Kerr nonlinearities, essentially ab- 
sorption free, that are up to 9 orders of magnitude larger 
than natural Kerr interactions (see Fig 1. (a)) [2(J. An- 
other distinctive feature of this scheme is that the non- 
linearity stays constant when several atoms interact with 
the same cavity mode. 

In this Letter we propose a new method for producing 
Kerr nonlinearities in cavity QED systems which (i) is 
experimentally less demanding, (ii) virtually absorption 
free, and (iii) produces nonlinearities comparable or even 
superior to the state-of-art EIT scheme [7[. By apply- 
ing suitable laser pulses at the beginning and end of the 
evolution of the proposed set-up, we obtain nonlinear in- 
teractions whose (iv) strength increases with increasing 
number of atoms interacting with the cavity mode, lead- 
ing to effective nonlinear interactions at least two orders 
of magnitude larger than previously considered possible. 
This brings closer to reality a number of proposals for 
quantum computation and communication based on pho- 
tonic nonlinearities. In particular, the estimated strength 
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FIG. 1: (a) EIT scheme considered in Refs. @, H, (b) dis- 
persive regime of the Jaynes-Cummings interaction, (c) large 
but noise nonlinearity scheme, and (d) scheme proposed in 
this paper. Coupling to the cavity mode is shown in green 
and to classical lasers in blue. 



of our nonlinearity is sufficiently large to allow the imple- 
mentation of quantum nondemolition measurements [llj , 
a photonic CNOT gate II 211 , and continuous variable en- 
tanglement distillation [131 ]. 

In our approach, the relevant atomic level structure, 
depicted in Fig. 1 (d), is a A system with two metastable 
states and an excited state. The cavity mode couples 
dispersively only to the — 2 transition and the levels 
and 1 are coupled via a far-detuned Raman transi- 
tion. Finally, the detuning associated with the lasers 
and with the cavity mode are assumed to be very dif- 
ferent from each other. This experimental setting can 
be implemented in several cavity QED systems, includ- 
ing atoms trapped in Fabry-Perot cavities 0] or opti- 
cal microcavities [3, EEI and quantum dots embedded 
in photonic band gap structures jl6| . In particular, our 
proposal is of considerable importance for the realization 
of many-body phenomena in systems consisting of ar- 
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rays of coupled cavities. There it was proposed 1(J 17 1 
that quantum-phase transitions for polaritons and pho- 
tons, as well as a photonic Mott-insulator [§] - in which 
photons would be frozen inside each cavity, mimicking a 
light crystal - could be experimentally observed provided 
strong photonic nonlincarities are available. 

To understand intuitively how our scheme works, let 
us start with the simplest system producing a nonlincar- 
ity: N two level atoms interacting dispersively with a 
cavity mode a (see Fig. 1 (b)). It is well known that 
if yfNg/A -c 1, where g is the Rabi frequency of the 
Jaynes Cummings interaction and A the detuning, then 
the system is described, including fourth order terms in 
the perturbation theory and neglecting non-energy pre- 
serving terms, by the following Hamiltonian: 



H eff = — aa ( S 



A 



Sn 



A 3 



a* aa' a(Soo — S11) 



where S aa :— \a-k)( a k\- Hence, if all the atoms are 

prepared in their ground states, we obtain a self Kerr 
nonlinearity of strength Ng 4 /A 3 . Despite the simplicity 
of this scheme, it has two major drawbacks. First, the 
obtained nonlinearities are at least one order of magni- 
tude smaller than in the EIT setting [20] . Second, due to 
the adiabatic condition s/Ng/A -C 1 characterizing the 
dispersive regime, the nonlinearity decreases as g/y/N. 

A simple solution to both these problems is to add a 
classical laser in resonance to the transition |0) — > |1) 
(see Fig. 1 (c)). Under the conditions \/Ng/A <C 1 
and \/Ng 2 /A <g; fi, the dynamics of the system is well 
described by 
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where S 3 := £f=i \+ k )(+ k \ - \- k )(- k \, with |±) = 
(|0) ± \1))/V2. Hence, preparing the atomic states in 
the superposition |— }, we find a nonlinearity as large as 
in the EIT setting [20] which does not decrease with the 
number of atoms. Unfortunately, all this improvement 
comes at the price of a prohibitively large increase in the 
losses of the system. Spontaneous emission will destroy 
the atomic superpositions, quickly decreasing the effec- 
tive nonlinearity and introducing noise to the system. 

In the scheme that we propose (see Fig. 1 (d)), we 
have a situation similar to the one discussed in the pre- 
vious paragraph. Indeed, the effective Hamiltonian will 
be identical to Eq. |T]). However, now the states |±) are 
metastable, which makes the system almost decoherence 
free. The full Hamiltonian of the system, in the interac- 
tion picture, reads: 
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where S 20 := £ fc |2 fc )(0 fc | and S 2 + := Ek\ 2 k)(+k\- As 
shown in Fig. 1 (d), (g, Ai) and (A, A 2 ) are the Rabi 
frequencies and detunings of the cavity-atom and laser- 
atom interactions, respectively. We are interested in the 
dispersive regime, characterized by: 
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Moreover, we assume that 



A^,VivA< |A 2 - Ai| 
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so that we can treat the processes driven by the cavity- 
atom and laser-atom interactions independently. Then, 
the excited state will hardly be populated, and we can 
adiabatically eliminate it, finding an effective Hamilto- 
nian for the two metastable states and the cavity mode. 
In turn, these will experience a.c. Stark shifts due to the 
interaction with the upper level. The effective Hamilto- 
nian, dropping out terms proportional to the identity, is 
given by 

g 2 9 
Hi = — a^aSoo + — (S w + S i), (5) 
Zai Z 

with G := 2A 2 /A 2 . If we now go to a second interaction 
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picture with respect to Hq — ^^a'a we find 

H int = _^_ flta {S+ _ + S _ +) + © 53; (6) 



where Sj 



S y = (S-\ ) t , and, as be- 



-*><- 

fore, S 3 := J2k=i l+fc)(+fe| _ |— fc)(— fc|- It follows that in 
the {|±)} basis the system can be viewed as an ensem- 
ble of two level atoms driven by a laser with a photon- 
numbcr-dependent Rabi frequency. If we consider the 
dispersive regime of this system, i.e. 
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the atoms prepared in the |— ) state will experience a 
Stark shift proportional to (a^a) 2 , which gives rise to 

the desired Kerr nonlinearity. The effective Hamiltonian 
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will be given by H e f f — —2- 
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i^na^aS^. Therefore, if we 
prepare all the atomic states in the |— ) state, one obtains 
an essentially absorption free Kerr nonlinearity [2]| given 
by 
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From Eq. ([8]) we can see that the strength of the effec- 
tive Kerr nonlinearity does not depend on the number 
of atoms, as we can tune Ai and independently. As 
such, the maximum achievable nonlinearity is limited by 
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g, which depends on the properties of the cavity and of 
the atoms employed. 

We now proceed to show that, in fact, our set-up can 
be modified to give nonlinear interactions which increase 

with N. The joint atomic operators ,5 |_, and S3 

satisfy the su(2) commutation relations: [S+-,S |_] = 

S3, [S3, S±zf:} = ±S± T . Defining the canonical trans- 
formation U = exp(na?a(S-\ — S [_)), with fi := 

g 2 /(Ai&) <C 1 [HI], we can use Hausdorff expansion 
{exp(xA)B exp(-xA) = B+x[A, B]+x 2 [A, [A,B}}/2+...) 
to obtain from Eq. © 

H rot := u\Hl nt )U^®S 3 -(^j i(ata) 2 S 3 

+ Q^) ^af(S + -+S- + ). (9) 

Suppose we had a way of generating H rot . If we pre- 
pared all the atoms in the |— ) state, the effective pho- 
tonic Hamiltonian would be given by the second term in 
Eq. © as long as 

This condition ensures that basically no transition from 
|— ) to |+) happens. Note that it is much less strin- 
gent than the one given by Eq. ([7]). In particular, set- 
ting 6 such that ( 5 2 /2Ai)/9 = iV- 1 ^ Eq. (33JJ) is 
satisfied for large N, and we obtain a nonlinearity of 
(y/Ng/A^N^^ga^aa^a. For instance, with N = 10 4 
and (y/Wg /Ai) = 0.1 , we obtain a nonlinearity equal 
to the Rabi frequency g, which is at least two orders of 
magnitude larger than in the single atom case. 

It is possible to realize the unitary operator V(t) = 
exp(itH rot ) for every t (see Fig. 2). We have that 
V(t) :— W exp(itH[ nt )U , hence it sufficies to show how 
to create the unitary U. Consider the Hamiltonian given 
by Eq. ([2]) when the classical lasers are swiched off, i.e. 
H = g(e- lAlt aS 20 + e iAl *a t 5 02 ). Suppose we apply to 
the atoms a fast laser pulse, described by the unitary 
operator M that generates the transformation 

|0fc)^(|0 fc )+i|l fe ))/V2 
|l fc )->(|0 fc )-i|l fc ))/V2, 

let Hamiltonian H run for a time t and apply the inverse 
transformation MT Then, the total evolution operator 
will be given by 

M*exp(ig f (e- iAlt 'aS 2 o + e lAlt ' a f S 02 )dt')M 
Jo 

= exp(i<? f V2g{e- lAlt ' a(S 20 + iS 2 %) + h.c.)dt'). 
Jo 

As \/Ng/A\ <C 1, the unitary operator above is well 
approximated by exp(i(<? 2 / A\)a) a(S+ — S 1 -)), which is 
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FIG. 2: Sequence of operations generating V(t). Blue and red 
areas correspond to the time the lasers with Rabi frequencies 
Q and A are on, respectively. The green area illustrates that 
the Jaynes-Cummings interaction between the cavity mode 
and the atoms is on at all times. The scales do not reproduce 
reality in general. For example, t will be much larger than 
I/O in most cases. 



U when t = I/O. The sequence of operations executing 
V(t) is shown in Fig. 2. 

In order to check the accuracy of this approach, we 
simulated it numerically for one and two atoms. We as- 
sumed that the laser pulse is generated by a 7r/2 pulse of 
the Hamiltonian H laser = £l(|0){lj + |1)(0|). We chosed 
fi = IOO5, which gives a t^/ 2 = 7r/(2fi). The other 
parameters are g = 10 8 s _1 , 50 = gN 1 ^ 3 , and Ai = 
lOyNg. With these values, the quality of the approxima- 
tion, which is determined by the value of the L.H.S. of Eq. 
(fTT)|) and y/Ng/Ai, should be constant. A good figure of 
merit in this respect is given by \(n\ ® (— |V(t)|— ) <8> |n)|, 
where V(t) is the total unitary operator formed by con- 
catenating the steps explained above (see Fig. 2 for a 
graphic description of the procedure), |n) is the n-th Fock 
state, and |— ) := | — 1) ® ... ® \~n) (see Fig. 3 (a)). In 
Fig. 3 (b) we plotted Re{(n\ <g> (-|V(t)|-) ® |n)) for 
(N,n) = (1,1), (1,2), (2,1), (2, 2), with the parameters 
g = 10 8 s" 1 , Ai = Wg, 6 = g, and Q = IOO3. We found 
a very good agreement with the dynamics of a pure pho- 
tonic nonlinearity of strenght k = N g 4 / 4AfQ , for which 
Re{(n\ (g) (-|V(*)|-) ® \n)) = cos{nn 2 t). 

Up to now we have discussed the case of self-Kerr non- 
linearities. In a cross-Kerr nonlinearity one optical field 
induces a Stark shift in another field proportional to the 
intensity of the former. In terms of two optical modes 
a and b, the Hamiltonian is given by H cross — va'aiwb. 
There are two possible generalizations of the scheme in- 
troduced here to produce large cross-Kerr nonlinearities. 

A shown in Fig. 4 (a), we may chose setting as in Eq. 
((21) , but now, we have a second cavity mode b coupled to 
the transition |1) — > |2) with Rabi frequency gb and De- 
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FIG. 3: (a) X = \{n\ ® {-\V(t)\-) ® |n)| versus t for (N,n) = 
(1,2) (blue) and (2,2) (green), with g = 10 V - 1), 59 = 
gN 1/3 , n = WOg, and A x = lOVTVg; (b) y = i?e((nj ® 
< — | V(t) I — > ® |n)) versus t for (JV,n) = (1, 1) (blue solid line), 
(1, 2) (blue dashed line), (2, 1) (red solid line), and (2, 2) (red 
dashed line), with g — 10 8 s _1 , Ai = lOg, = g, and Q = 
100 5 . 



(a) 



(b) 





FIG. 4: Level structures for generating cross-Kerr nonlinear- 
ities. 



tuning Ai b . As shown in Ref. [19( , this could be achieved 
using two degenerate cavity modes with orthogonal po- 
larization. Carrying over the same analysis we did for 
the one mode case, we can find in this case 
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Hence, in addition to a self-Kerr nonlinearity for each 
mode, we find a cross-Kerr nonlinearity between them. 

In the second case, we consider the level structure 
shown in Fig. 4 (b) . We have two modes a and b coupled 
to the transition |0) — * |2) with detunings Ai ± 8, re- 
spectively. As discussed in Ref. [14j , this set-up could be 
realized in a toroidal microcavity, where a and b are the 
normal modes of the clock and counter-clockwise prop- 



agating modes, and 8 is the rate of tunneling between 
those two. In this case, the effective Hamiltonian in this 
case can be found to be 
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Summarizing, we have presented a new proposal for 
producing giant Kerr nonlinearities in cavity QED sys- 
tems. Our approach could be used to generate nonlin- 
earities at least two orders of magnitude larger than pre- 
viously considered possible. These, in turn, could be ap- 
plied to the implementation of several quantum infoma- 
tion processing tasks with photons, as well as paving a 
way to the observation of quantum many-body phenom- 
ena in coupled array of cavities. 
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